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We investigate the prospects of indirect and direct dark matter searches within the minimal su¬ 
persymmetric standard model with nine parameters (MSSM-9). These nine parameters include 
three gaugino masses, Higgs, slepton and squark masses, all treated independently. We perform a 
Bayesian Monte Carlo scan of the parameter space taking into consideration all available particle 
physics constraints such as the Higgs mass of 126 GeV, upper limits on the scattering cross-section 
from direct-detection experiments, and assuming that the MSSM-9 provides all the dark matter 
abundance through thermal freeze-out mechanism. Within this framework we find two most prob¬ 
able regions for dark matter: 1-TeV higgsino-like and 3-TeV wino-like neutralinos. We discuss 
prospects for future indirect (in particular the Cherenkov Telescope Array, CTA) and direct detec¬ 
tion experiments. We find that for slightly contracted dark matter profiles in our Galaxy, which 
can be caused by the effects of baryonic infall in the Galactic center, CTA will be able to probe a 
large fraction of the remaining allowed region in synergy with future direct detection experiments 
like XENON-1T. 


I. INTRODUCTION 

Identifying the particle nature of dark matter is one 
of the most pressing goals of modern astrophysics and 
cosmology. If dark matter is made of weakly interacting 
massive particles (WIMPs) 2 [2], where the relic dark 
matter abundance is naturally explained with the ther¬ 
mal freeze-out mechanism, a particular promising avenue 
is the search for signatures from the self-annihilation of 
dark matter particles in gamma rays [3] . The required av¬ 
erage velocity-weighted annihilation cross-section during 
freeze-out is of the order of (av) = 3x 10 -26 cm 3 s _1 [TJ 
|4j. Upper limits from modern gamma-ray instruments, 
such as the Large Area Telescope (LAT) aboard Fermi 
satellite, start to exclude this canonical annihilation 
cross-section for WIMP masses below 100 GeV [3- At 
the same time, the hunt for signatures of physics beyond 
the standard model of particle physics, especially of su¬ 
persymmetry, is ongoing at the Large Hadron Collider 
(LHC). Supersymmetric models provide interesting can¬ 
didate particles for WIMP dark matter (in most cases the 
lightest neutralino). Although there has been no claim of 
positive signatures of supersymmetry yet, some outstand¬ 
ing conclusions can be drawn from several measurements 
at the LHC. One of them comes from the mass measure¬ 
ment of the Higgs boson of about 126 GeV [3IZj. 

Since particle masses and couplings are subject to 
loop corrections, the Higgs mass measurement can pro¬ 
vide indirect information about yet undiscovered parti¬ 
cles. Therefore, this measurement already significantly 
constrains the parameter space of supersymmetric mod¬ 
els. Reference [8] studied the implications of the Higgs 
mass measurement for the minimal supersymmetric stan¬ 
dard model (MSSM) with five parameters (the con¬ 


strained MSSM), and found that the posterior proba¬ 
bilities for these parameters are narrowly distributed. 
Reference [9] further extended the model to include the 
non-universality in gaugino and Higgs masses, and found 
qualitatively similar results. The most probable masses 
of neutralino dark matter were found to be around 1 TeV, 
for higgsino dark matter, and 3 TeV, for wino dark mat¬ 
ter. 

The wino dark matter case is already in some ten¬ 
sion with searches for gamma-ray lines [T0] - [T2] . The 
wino dark matter annihilation cross-section is signifi¬ 
cantly larger than the canonical value for WIMPs due 
to a non-perturbative effect known as Sommerfeld En¬ 
hancement (SE) flBlHTj . The SE calculation for heavy 
WIMPs is subject to large logarithmic corrections due to 
the large hierarchy between the dark matter mass and 
the W-boson mass. Reference m showed that full one- 
loop computation makes the cross-section smaller up to 
about 30% with respect to the SE correction at tree level. 
Even more precise computations using Soft Collinear Ef¬ 
fective Theory were recently presented (Tf)U221 . Refer¬ 
ence [22] found that when calculating leading-log semi- 
inclusive rates, the effect of higher order corrections are 
very modest. 

There are several studies in the literature that have in¬ 
corporated the impact of the Higgs mass measurement, 
and of XENON-100 [23] and LUX (23] bounds, on the 
WIMP-nucleon spin-independent scattering cross-section 
on constrained MSSM scenarios (see, e.g., Refs. [8] 1011251- 
1501). Various statistical approaches have been used to 
infer the most probable regions of these scenarios. The 
parameter space is often restricted to energies below a 
few TeV, according to what is expected from “natural” 
supersymmetry. Interestingly, when performing a proper 
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Bayesian analysis, the fine-tuning penalization arises au¬ 
tomatically from very basic statistical arguments (the 
Bayesian version of “Occam’s razor”), allowing to ex¬ 
plore larger regions of the parameter space while taking 
the notion of naturalness automatically into account (see, 
e.g., Ref. [131]). 

Reference [5] studied the non-universal higgsino and 
gaugino masses model within the Bayesian framework. 
The authors showed that the most probable regions for 
the neutralino dark matter are around 1 TeV (higgsino 
dark matter) and 3 TeV (wino dark matter), the high 
masses being mostly due to the Higgs mass and the relic 
density constraint. Although this leaves a very large 
portion of the favoured supersymmetric parameter space 
outside the reach of LHC, prospects for future dark mat¬ 
ter experiments were shown to be very promising, in par¬ 
ticular for XENON-1T [32] for direct detection, and for 
the Cherenkov Telescope Array (CTA) [33] for indirect 
searches. Recently, Ref. [33] presented a first study of the 
detection prospects for CTA and XENON-1T for neu¬ 
tralino dark matter in the 19-parameter phenomenologi¬ 
cal MSSM. In this study, the SE was effectively included 
as an extrapolation of the results from Refs. Ha mi. 

In the present paper, we study consequences of MSSM 
models with extended nine parameters (MSSM-9), and 
prospects for indirect and direct dark matter searches. 
In addition to the non-universality of the gaugino and 
Higgs masses studied in Ref. [9], we investigate the non¬ 
universality of masses and the trilinear couplings in the 
sfermion sector (i.e., they are independent between slep- 
tons and squarks). This treatment is more general be¬ 
cause constraints from the LHC mainly affect the squark 
and gluino sector and do not directly reflect on the 
sleptons. Since the preferred dark matter masses are 
at 1 TeV and above, we discuss constraints from the 
Fermi gamma-ray satellite and current generation of 
Cherenkov telescopes, in particular HESS [35] , As men¬ 
tioned above, the wino dark matter around 3 TeV is 
subject to the SE correction of the annihilation cross- 
section m mj m on, and, therefore, we calculate the 
SE point-by-point in the scan. In addition, the annihi¬ 
lation of wino dark matter also yields strong gamma-ray 
line signals, which are tightly constrained by the HESS 
observations of the Galactic center [35]. Note that our 
results are based on a full numerical study of MSSM pa¬ 
rameter scan, and does not rely on a simplifying assump¬ 
tion of pure wino case as adopted in Refs. |10] - fT2] . We 
will discuss prospects for future indirect and direct detec¬ 
tion experiments, most notably for CTA that will have 
an excellent sensitivity for gamma rays above 100 GeV 
[37H!Ij. 

The paper is organized as follows. In section II, we 
discuss the model and the adopted scanning technique. 
We will summarize the results of the scan in section III, 
and show in section IV prospects for future indirect and 
direct detection experiments. In section V, we give our 
conclusions. 


II. MSSM MODELS AND HIGGS MASS 


The discovery of the Higgs boson [313 has completed 
the picture of the standard model of particle physics, 
which has proven an extremely good description of par¬ 
ticle physics up to the TeV scale. Beyond the crucial 
importance of this discovery by itself, this result has far- 
reaching consequences for well-motivated candidates of 
physics beyond the standard model, such as supersym¬ 
metry, and in particular for the MSSM. 

The rather high reported Higgs mass m^ shifts the 
scale of supersymmetry to higher values. In the MSSM, 
the tree-level Higgs mass is bounded by the mass of the 
Z-boson, and, therefore, large radiative corrections are 
needed in order to reconcile theory and experiment. An 
approximate analytic formula for m^ [33 @3] reads 
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where tan /3 is the ratio between the vacuum expecta¬ 
tion values of the two Higgs doublets v u = and 

Vd = (H < d)i v 2 = v 2 + m t is the top running mass, 
and Msusy represents a certain average of stop masses. 
X t = A t — p cot (3, where /i is the Higgs mass term in 
the superpotential, and A t is the trilinear stop coupling, 
both at the electroweak breaking scale. The first term 
of Eq. ([!]) is the tree-level Higgs mass, while the sec¬ 
ond two terms are the dominant radiative and thresh¬ 
old corrections. Note that the radiative corrections grow 
logarithmically with the stop masses while the threshold 
corrections has a maximum for X t = ±\/6Msusy- To 
achieve m/, ~ 126 GeV, one typically needs stop masses 
larger than ~3 TeV, unless X t is close to its maximum 
value. 

In order to evaluate the sensitivity on dark matter in 
a more generic context, we parameterize the MSSM with 
10 fundamental parameters at the gauge coupling uni¬ 
fication scale. After requiring the correct electroweak 
symmetry breaking, we end up with 9 effective parame¬ 
ters: 

js, Mi, M 2 , M 3 , mg, m H , Ag, Ag, tan/3, sgn(/it)} (2) 


where s represents the SM nuisance parameters, 
Mi,M 2 , M 3 are the gaugino masses, ?ng, mg, mH are the 
squark, slepton, and Higgs masses (mn = = m H d ), 

and Ag and Ag are the squark and slepton trilinear cou¬ 
plings. The sign of p is fixed to +1. All the soft param¬ 
eters defined at gauge coupling unification scale, except 
for tan/3 and the SM nuisance parameters. Compared 
with Ref. (9], we further generalize the sfermion sector, 
by adopting independent values for squarks and sleptons. 

We perform a Bayesian analysis to generate a map of 
the relative probability of different regions of the param¬ 
eter space. In doing so, the global likelihood is defined as 
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Observable 

Mean value Uncertainties 

fi a (exper.) r (theor.) 

Ref. 

M w [GeV] 

80.399 

0.023 

0.015 

nr 

sin 2 6 e f f 

0.23153 

0.00016 

0.00015 

m 

BR (B ->• X s7 ) x 10 4 

3.55 

0.26 

0.30 

m 


1.04 

0.11 

- 

m 

BR(S„->riO SM 

A 0 _ x 10 2 

1.63 

0.54 

- 

ESI 

3.1 

2.3 

- 

m 

BR (B^Drv) w , n 2 

BR (B^-Deu) A 

41.6 

12.8 

3.5 

ESI 

Rl23 

0.999 

0.007 

- 

m 

BR(_D S —>■ tv ) x 10 2 

5.38 

0.32 

0.2 

:45> 

BR(D S —> fj,v) x 10 3 

5.81 

0.43 

0.2 

HU 

BR(D —> /iv) x 10 4 

3.82 

0.33 

0.2 

m 

Q. x h 2 

0.1196 

0.0031 

0.012 

eh 

rrih [GeV] 

125.66 

0.41 

2.0 

eh 

BR (B s —> 

3.2 x 10 

" 9 1.5 x 10~ 9 

10% 

EH 


Limit (95% CL ) 

Ref. 

Sparticle masses 

ST 

m x - a xN 

As in Table 4 of Ref. ]53j. 

XENON100 2012 limits 

m 


TABLE I. Observables used for the computation of the like¬ 
lihood function 

a multiplication of individual likelihood functions, where 
for each quantity we use a gaussian function with mean 
/i and standard deviation s = \Jo' 1 ± t 2 , where a is the 
experimental uncertainty and r represents our estimate 
of the theoretical uncertainty. For upper and lower limits 
we use a gaussian fuction to model the drop in the like¬ 
lihood below or above de experimental bound. The ex¬ 
plicit form of the likelihood function is given in ref. |53| . 
including in particular a smearing out of experimental 
errors and limits to include an appropriate theoretical 
uncertainty in the observables. We take into consider¬ 
ation all the available particle physics data described in 
table[TJ including the Z mass, which is effectively included 
adding a Jacobian factoiQ electroweak precision mea¬ 
surements [53] , B-physics observables [35] [38] [33] [53] [55] , 
the Higgs mass mm, and constraints on the WIMP- 
nucleon scattering cross-section by XENON-100 [33] and 
LUX [33] . In addition, we assume a scenario with a sin¬ 
gle dark matter component that is produced thermally in 
the early Universe, by including the measured relic den¬ 
sity according to the Planck results ;5BJ. For the relic 
density and (av) computation, we take the SE into ac¬ 
count by creating a grid of the enhancement in the M 2-/1 
plane using the Hryczuk et al. computation method im¬ 
plemented in DarkSE ca ng. For the computation of 
(av) in the present day we implemented a function in 
DarkSE to extract the enhancement for v = 1CU 3 from 
the Hryczuk et al. computation (we validated the results 
with the pure-wino case showed in [18]). 


1 The Jacobian factor arise after perform a change of variable 

{yi , //. B } —» {mi, M z , tan /?} and integrating M z in the posterior 

probability density function, as explained in detail if Ref. ED. 


Parameters 

Range scanned 

Mi, M 2 , M 3 [GeV] (-10°, 10°) 

[GeV] (10, 10 6 ) 

Al,A l Q [GeV] 

(-10 6 , 10 6 ) 

tan/3 

(2, 60) 

sg n (t0 

+1 


TABLE II. Ranges of model parameters adopted in the scan. 



Gaussian prior Range scanned 

ref. 

M t [GeV] 

173.2 ±0.9 (167.0,178.2) 

EH 

rrib(m b ) MS 

[GeV] 4.20 ± 0.07 (3.92, 4.48) 

m 

[a ern (M z ) MS ]~ 1 127.955 ± 0.030 (127.835, 128.075) [57] 

a s (M z ) MS 

0.1176 ±0.0020 (0.1096, 0.1256) 

EH] 


TABLE III. Nuisance parameters adopted in the scan. 

For the priors of the parameters, we adopted both stan¬ 
dard and ‘improved’ log priors (S-log and I-log, respec¬ 
tively), defined in Refs. [8, 9j, with ranges described in 
table [TT] In the case of the I-log priors, we effectively 
assume that the parameters are associated with a com¬ 
mon scale (motivated by a common underlying super- 
symmetry breaking mechanism). In the following, we 
show results for the I-log priors only, but note that those 
obtained with S-log priors are very similar, since the 
data, in particular the relic density constraint and the 
Higgs mass, turns out to be very constraining. For the 
SM parameters, we used gaussian priors described in ta¬ 
ble |III| For the numerical analysis, we use the Super- 
BayeS code [59], which uses the nested sampling algo¬ 
rithm implemented in Multinest [801 , and integrates Soft- 
Susy [63], SusyBSG [53], Superlso [533, DarkSusy 53] . 
MicrOMEGAs [85], and DarkSE [TT] for the computa¬ 
tion of the experimental observable. The full likelihood 
function is the product of the individual gaussian likeli¬ 
hoods associated to each piece of experimental data. In 
particular, for XenonlOO we use the likelihood defined in 
[66] . The LUX limit is applied as a step function only for 
the 2a confidence level points. For a more detailed ex¬ 
planation of the Bayesian analysis relevant for the results 
in this paper, we refer the reader to Ref. 9]. 

III. RESULTS OF THE SCAN 

In Fig. [I] we show two-dimensional contours that rep¬ 
resent 68% and 95% credible regions of the most rele¬ 
vant parameters for CTA: the dark matter neutralino 
mass and annihilation cross-section. The posterior has 
two peaks in the mass distribution. The largest peak lo¬ 
cates around 1 TeV, where the neutralino mostly con¬ 
sists of higgsino. There is a weaker peak around 3 
TeV, where it is mostly wino. The wino dark mat¬ 
ter features significantly larger annihilation cross-section 
around 1CV 24 cm 3 s _1 due to SE correction. Since the 
SE is a non-relativistic effect causing the distortion of the 
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FIG. 1. The contours represent 68% and 95% posterior proba¬ 
bility credible regions. Colored points reproduce all the exper¬ 
imental observables within 2<r of conhdence level. The cyan 
diamond represent the pure higgsino case from inn, and the 
blue triangle the pure wino case from 1101 118} . The color¬ 
coding indicates the branching fraction into W + W~ final 
states. The green lines show the sensitivity of CTA as de¬ 
rived in Ref. [40] for bb final states and lOOli observation time 
(assuming 1% systematics). The purple lines show the HESS 
Galactic center limits for generic hadronic final states [68] . 
adopting an Einasto (solid), a contracted NFW (dotted), and 
a shallow NFW (dashed) profiles The brown dotted line shows 
the Fermi-LAT limit from the analysis of the dwarf spheroidal 
galaxies [69j . The CTA and Fermi-LAT limits correspond to 
XiXi —t bb, but limits for W + W~ final states are very similar. 


wave functions, it is more efficient for heavier particles. 
We note that these two regions correspond to those found 
in Ref. |9] with a seven-parameter MSSM study. In fact, 
the most probable regions in the posterior distributions 
for the mass of the lightest particles (the neutralino) are 
only mildly changed compared to [0]. This shows the 
robustness of the procedure against the number of pa¬ 
rameters. 

Figure [T] also shows, as points, regions in the param¬ 
eter space that reproduce all experimental observables 
within 2cr of confidence level. We remind that the pos¬ 
terior probability distribution function (PDF) shows rel¬ 
ative probabilities within a model, given the experimen¬ 
tal data, under the hypothesis that the model is cor¬ 
rect. The 68 % and 95% credibility regions show that it 
is much more likely to find neutralinos with a mass of 
~1 TeV and ~3 TeV, however, these contours not neces¬ 
sarily cover all the regions that respect the experimental 
observables. The scattered points outside the contours 
show regions that require more tuning to reproduce the 
experimental observables and, therefore, their integrated 
probability is small. These less probable regions, with 


FIG. 2. The same as Fig. [T] but for annihilation into 
monochromatic photons. The cyan diamond represent the 
pure higgsino case from ,10, ;T8], and the blue triangle the 
pure wino case na. Colors indicate the higgsino fraction of 
the lightest neutralino. The green lines show the line sensi¬ 
tivity of CTA as derived from ED], while purple lines are the 
HESS limits |7l], all for the Galactic center. 


dark matter between 1 TeV and 3 TeV, correspond to 
wino-higgsino and wino-bino neutralinos. 

An additional region around hundreds GeV, corre¬ 
sponding to bino-like neutralino, is not show in the figure. 
Unlike higgsinos and winos, bino neutralinos can not self- 
annihilate, therefore, a specific mass relation with other 
mass eigenstate is required to have an efficient enough 
annihilation to reproduce the correct relic density. For 
example, a bino quasi-degenerate with the stau, or a bino 
mass equal to half of the lightest Higgs or pseudo-scalar 
mass. On the other hand, as we mention above, unless 
we are in the maximal mixing scenario (which is also sub¬ 
ject to certain tuning, see Ref. (72]) the Higgs mass mea¬ 
surement tend to push the spectrum to higher masses. 
For few hundreds GeV neutralinos, a fine-tuning is nec¬ 
essary to reproduce the Higgs mass and the relic density. 
Therefore, this region has small statistical weight and is 
not well explored in our scan. For this reason we show 
m x o larger that 500 GeV in our figures. In any case, few 
hundreds GeV neutralinos will be potentially tested by 
the LHC. 

The annihilation of wino-like dark matter produces rel¬ 
atively strong line signatures. Given that the total anni¬ 
hilation cross-section is very large, as shown in Fig. [T] the 
resulting annihilation cross-section into 77 or 7 Z chan¬ 
nel is correspondingly large. Figure [ 2 ] shows the contours 
of mass and (( 7 v ) 77 + {av) 1 z/‘2, which is the relevant 
quantity for line searches with Cherenkov telescopes since 
the energy splitting between these two modes are smaller 
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FIG. 3. Similar to Fig. [I] but showing the total annihila¬ 
tion cross-section against the WIMP-proton spin-independent 
cross-section. The blue triangle and the cyan arrow are the 
theoretical values of the pure wino and pure higgsino case 
from Refs. [73ll74j . The arrow indicates that the value of rjp 1 
is a theoretical upper limit. The color-coding indicates the 
higgsino fraction. Green lines show the sensitivity of CTA 
for Einasto (solid) and contracted NFW (dotted) profiles for 
the Galactic center [40]. Cyan and pink lines represent LUX 
limit [24] and XENON-1T sensitivity [32]. Both the limit and 
sensitivity lines assume m x o = 1 TeV. CTA sensitivity lines 

correspond to XiXi ~ 1 ► bb. 


than their typical energy resolution (e.g., Ref. HD]). 

The colored points in Fig. [l] show the branching frac¬ 
tion of the lightest neutralino annihilating into W + W~. 
For the higgsino-like neutralino (m x o ~ 1 TeV), the an¬ 
nihilation to ZZ becomes very important, while for the 
wino-like neutralino (m x o ~ 3 TeV), the annihilation to 
W+W~ is the dominant channel. The colored points 
in Fig. [2] shows the higgsino fraction of the lightest neu¬ 
tralino. Note that points inside the 95% contour are dom¬ 
inantly higgsino-like and wino-like neutralinos. However, 
by comparing those points with the pure higgsino (cyan 
diamond) and pure wino (blue triangle) cases, it is clear 
that the 95% contour encloses points with small mixing 
that could have different characteristics. In particular, 
in the wino-like region with (av) 11 J rl z /2 > 10“ 27 cm 3 /s 
and to. 0 ~ 2.5 TeV, the co-annihilation with sfermions 
plays a very important role in setting the correct relic 
density. 


IV. GAMMA-RAY UPPER LIMITS AND 
PROSPECTS 


The Galactic center is one of the most promising places 
to search for signals from WIMP annihilation, as it is 


locally the densest region of dark matter, even though it 
is often challenging to address astrophysical foregrounds 
(see, e.g., Ref. [75]). In general, the gamma-ray intensity 
from neutralino annihilation towards a direction 1 p away 
from the Galactic center is given by 


Ij (4 ; -®7 ) 


(fTu) dN ltSi nn 


8nm 2 0 

Xi 


dEry 


J(4 0 ) 


J(4’) = 


dlp 2 x 0 (r[l,ip ]), 


(3) 


where dN 1:ann /dE^ represents the annihilation spec¬ 
trum, r{l,ip) 2 = Tq + l 2 — 2r 0 l cosip, tq = 8.5 kpc is 
the Galactocentric radius of the solar system, and Z max 
corresponds to the virial radius of the Milky Way. The 
dark matter density profile in the Galactic halo p x .0 is 
widely assumed to be given by phenomenological fits to 
results of dark-matter-only N-body simulations, such as 
the Navarro-Frenk-White (NFW) [7Bj or Einasto [771178] 
functions. We note, however, that these profiles are ob- 
servationally confirmed only at relatively large radii, and 
far less constrained in the inner regions tzE where the 
strongest annihilation signals are expected. We will com¬ 
ment again on this point below. 

a. Current Status. Figures [l] and [5] show the pre¬ 
dicted annihilation cross-section into continuum pho¬ 
tons (dominated by W + W~, ZZ and qq final states) 
and gamma-ray lines, respectively, compared to differ¬ 
ent experimental limits and reaches. Our fiducial density 
profile is given by an Einasto profile (with parameters 
a = 0.17, pQ = 0.4GeV/cm 3 and r s = 20kpc). The 
current upper limits on the Galactic center from HESS 
searches for gamma rays from bb final states (Fig. 0 EH 
ED and for gamma-ray lines (Fig. ©in] are already very 
tight. We find that the wino dark matter region around 
3 TeV is almost completely excluded by the HESS up¬ 
per limits. This is in agreement with the findings of 
Refs. [Tdrl 12 1 . but interpolating in the av enhancement 
from a M 2 - p grid instead of extrapolated from exist¬ 

ing pure wino calculations as is in Ref. [34] . 

We note, however, that upper limits are still subject 
to uncertainties mainly related to the density profile m- 
While state-of-the-art N-body simulations prefer either 
NFW or Einasto profiles, baryonic effects can potentially 
modify the density profiles significantly. For example, 
baryonic adiabatic contraction can compress the dark 
matter profiles, while supernova feedback can cause the 
profiles to be shallower (see, e.g., Refs. j8UH85| ). To illus¬ 
trate this point, we show in Figs. [T] and [ 2 ] how the HESS 
limits weaken when a shallower dark matter profile is 
adopted. To this end, we use a generalized NFW pro¬ 
file with an inner slope of 7 = 0.7 (and r s = 20 kpc, 
Pq = 0.4GeV/cm 3 ), which is still in agreement with 
kinematic and microlensing observations [79]. In this 
case, the limits indeed are weakened and part of the wino 
best-fit region is still allowed. A similar effect will occur 
for cored profiles. 
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In Fig. |TJ we also show the Fermi -LAT limits from the 
observation of dwarf spheroidal galaxies from Ref. :69j, 
which already include the uncertainties in the dark mat¬ 
ter profile and can be hence considered robust (i.e., this 
represent the upper end of the uncertainty band). They 
exclude most of the wino parameter space. 

b. Prospects for CTA. Current constraints leave the 
1-TeV higgsino dark matter as most interesting dark mat¬ 
ter candidate in the MSSM-9. The CTA sensitivities for 
both the total annihilation cross sectionj40] and for the 
gamma-ray lines EH! are shown in Figs, [l] and [2j respec¬ 
tively. These figures show that, for standard Einasto pro¬ 
files, it will be challenging for CTA to reach the 1-TeV 
higgsino parameter space, unless background systematics 
are under control at the sub-percent level [JO]- However, 
as explained above, baryonic effects could potentially in¬ 
crease the chances for a CTA discovery of higgsino dark 
matter as baryons can drag dark matter towards the 
Galactic center during their cooling, leading to a more 
cuspy profile [50]. To illustrate this effect, we addition¬ 
ally show in Fig. ljthe reach of CTA when a slightly con¬ 
tracted NFW profile, with an inner slope of 7 = 1.3 (and 
otherwise parameters as above), is adopted. In this case, 
CTA has the potential to rule out (or discover) a large 
part of the best-fit higgsino dark matter region. Our re¬ 
sults are less stringent than the results found in Ref. , 
since our estimates for the CTA sensitivity also include 
an estimated 1% systematic uncertainty. Note also that 
a more accurate treatment of the cosmic ray background 
in Ref. [JT] leads, for dark matter masses above 1 TeV, 
to slightly less stringent projected limits than what we 
show here, within a factor of two. 

c. Direct Detection. Fig. [3] shows the most proba¬ 

ble regions plotted for the annihilation cross-section and 
the spin-independent scattering cross-section tjp 1 at tree 
level. To give a reference of the size of for pure 
cases, we also show higgsino neutralino (cyan diamond) 
and wino neutralino (blue triangle) one-loop computation 
performed by In the pure higgsino case, the 

perturbative QCD and hadronic input lcr-uncertainties 
allow only to set a maximum value for a ^ 1 , which is rep¬ 
resented in the figure by an arrow (^] Hence, the mixing 
of the neutralino has a crucial role for the computation 
of (jp 1 . The colored points in Fig. ^represent the hig¬ 
gsino composition of the lightest neutralino, and show 


2 The tree level scattering of the neutralino with the nucleon, 
through the Higgs boson, requires a neutralino with non- 
negligible higgsino-bino or higgsino-wino mixing. The other pos¬ 
sibility is the scattering via squarks, in that case squarks should 
be light enough to give a sizable contribution. 

3 For the computation of a SI we have used the s-quarks nucleon 

form factor derived from measurements of the pion-nucleon sigma 

term |84| . On the other hand, refs. EH ED use the lattice calcu¬ 

lation value. Using the lattice value the tree level cross section 

drops by a factor of 4 ref El- 


how the value of a^ 1 decreases with the higgsino frac¬ 
tion in the wino-like region. We remind the reader that 
this computation was performed at tree level in our scan. 
However, almost pure wino points appear in the region 
where cr SI is smaller that ~ 10 -11 pb and av larger than 
10~ 25 cm 2 3 * * * * s _1 in Fig. [3| Comparing these a SI values 
with the one of the pure wino case (blue triangle), it 
is clear that the one-loop contribution is the dominant 
one. Therefore, we would expect that after including 
higher order corrections those points will get a SI values 
of ~ 10” 11 pb. 

Figure[3]also shows the sensitivities of XENON-1T [32] 
and CTA [JO] (both for m % o = 1 TeV), showing that 
both direct and indirect searches are very important for 
the potential discovery of TeV dark matter, that is, at 
the moment, the most probable solution in the context 
of MSSM-9. Note also that the region around the almost- 
pure and pure wino and higgsino neutralinos will be 
probed by CTA only. 


V. CONCLUSIONS 

We studied the prospects for indirect and direct dark 
matter searches in context of the MSSM-9 by means of 
a Bayesian Monte Carlo scan. We find as the two most 
likely regions the 1-TeV higgsinos and 3-TeV winos dark 
matter. Current limits from dwarf spheroidal observa¬ 
tions as well as observations of the Galactic center with 
Fermi -LAT and HESS exclude almost all the models with 
wino-like dark matter, even for flattened profiles of the 
dark matter halo. However, models with 1-TeV higgsino- 
like dark matter remain unconstrained. We find that 
for regular dark matter profiles, it will be challenging 
for CTA to probe the higgsino dark matter parameter 
space, both for continuum and gamma-ray line searches. 
However, a mildly contracted profile would improve the 
prospects significantly, and make most of the higgsino 
dark matter parameter space testable in the upcoming 
years, providing complementary constraints to future di¬ 
rect detection experiments like XENON-1T. 
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